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We calculate the electric flux and magnetic monopole current distribution in the 
presence of a static quark-antiquark pair for SU(2) lattice gauge theory in the max- 
imal Abelian gauge. The current distribution confines the flux in a dual Abrikosov 
vortex whose core size is comparable to the flux penetration depth. The observed 
structure is described by a dual Ginzburg-Landau model. 
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I. INTRODUCTION 



One explanation for the absence of free quarks is that the QCD vacuum naturally expels 
color-electric flux in a manner similar to the expulsion of magnetic flux in a superconductor. 
In this scenario the vacuum must contain objects that respond to a color-electric field 
by generating screening currents to confine the electric flux in a narrow tube similar to 
the Abrikosov flux tube produced in a superconductor || . While this can be achieved in 
Abelian models [[!]] by adding additional Higgs fields, in the dual superconductor model |2| 
dynamically generated topological excitations (magnetic monopoles) generate the currents. 
In 4D U(l) lattice gauge theory there is considerable indirect evidence from "bulk properties" 
of the vacuum such as the monopole density f|, monopole susceptibility ||, and static 
quark potential that Dirac magnetic monopoles are associated with the phenomenon of 
confinement. 

In a recent paper JFj we reported the first direct evidence that the monopoles in U(l) 
lattice gauge theory actually react to the electric flux from a static quark-antiquark pair 
by producing a solenoidal current distribution to confine the flux to a narrow tube. In the 
confined phase the curl of the monopole current and the local electric field obey a dual version 
of the London equations || for a superconductor. We also found evidence for quantization of 
the electric fluxoid ||. The flux tube that formed between a static qq pair was structureless 
with no "normal core" visible on the scale of a lattice spacing. In the deconfined phase none 
of this behavior was observed. 

While the origin of confinement in U(l) lattice gauge theory is fairly clear, progress in 
understanding confinement in non-Abelian SU(N) theories has been slow. One promising 
approach || is to fix the non-Abelian degrees of freedom in the maximal Abelian gauge P1,|10 



leaving a residual U^l)^ -1 gauge freedom, with (iV-1) species of U(l) Dirac monopoles. The 
monopoles have been observed to be [p|Hj3|| abundant in the confined phase and dilute in 
the (finite temperature) unconfined phase. We present here the first direct evidence that 
a dual Abrikosov vortex also forms in SU(2) lattice gauge theory with static quarks. We 
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investigate the structure of the Abrikosov vortex and find that the flux tube in the present 
case has a normal core of size comparable to the flux penetration depth. Our results are 
consistent with a dual version of the Ginzburg-Landau model of superconductivity j3|. A 
preliminary report of our work was presented at the LATTICE 92 conference ||14|| . 



II. SIMULATIONS 

Our simulations were performed on a 13 3 x 14 lattice with skew-periodic boundary con- 
ditions. Each link from point s in the /i direction carried an SU(2) element U(s, fi) and the 
plaquette operators P^ u {. s ) were formed in the usual fashion as a directed product of link 
variables. We used a standard SU(2) Wilson action 



S = ^(l--ReTr P^a)) . (1) 
We generated configurations distributed according to exp(— S) using a combination of Monte 



Carlo [jHJ and overrelaxation [16|. Simulations were performed for (3=2 A and 2.5 with an 
initial 1000 sweeps to thermalize followed by 50 sweeps between measurements. 

We converted our configurations to the maximal Abelian gauge || by finding the gauge 
transformation that maximized the quantity |9| 

R = J2 Tr [cx 3 U(s,fi)a 3 U\s,fi)}, (2) 

S,fl 

which is equivalent to diagonalizing 

x ( s ) =^2[ u ( s > ^)^U ] (s, fi) + U(s - n, /j)a 3 U j (s - //, n) (3) 

at each site s. To measure the gauge fixing in the simulations we used the lattice sum of the 
magnitude of the offdiagonal component of X(s), \Z\ 2 = J2 S |A(s)i2| 2 . Typically we needed 
about 650 gauge fixing sweeps to attain \Z\ 2 w 10 _5 /site for (3 = 2.4. Three gauge fixing 
methods were used: (1) generating and accepting random local changes only if R increased, 
(2) locally maximizing R exactly at alternate sites, (3) applying overrelaxation using the 
square of the gauge transformation of method (2) to sample configurations better. 



III. MEASUREMENTS 



After gauge fixing, the Abelian U(l) link variable is given || by the phase of the diagonal 
component of the SU(2) link variable, u(s,fi) = U(s, f/,)u/\U(s, //)n|. We construct the 
Abelian plaquette variables p^ v {s) from the u(s, /i) in the usual fashion. The one important 
difference from our U(l) study is that we have no local U(l) action here, so we cannot 
identify the magnitude of the Abelian charge e. In the U(l) case we divided the plaquette 
by e in order to get the electric field E, and we multiplied the integer valued lattice operator 
for V x J M by cm = 2ir/e to normalize the magnetic current. In the present case we 
determine instead the electric "force" eE and the magnetic number current 2it Jm / 'e-M — &Jm 
as variables here, so both quantities are known only up to an overall common factor. We 
shall continue to refer to them as the electric field and magnetic current in keeping with our 
previous analysis. 

Since the electric field and monopole current operators are vectors, they will average 
to zero unless they are correlated with a Wilson loop W representing the current from a 
qq pair; a 3 x 3 loop in the z-t plane was used in the simulations. Therefore, averages of 
observables are computed as (6) = Tr{exp(— S)W Q} /Tr{exp(— S)W}. 

The electric field operator is given by a 2 8^ = Im p M4 , where a is the lattice spacing. With 
our choice for the orientation of the Wilson loop only the z-component of the electric field 
has a nonzero average. In Fig. |l](a) we show the operator for the electric field £ z , given by 
a z-t plaquette, bold line for fixed time. 

The magnetic monopoles are identified using the DeGrand-Toussaint 0] construction. It 
is convenient to associate the monopole current density in each three-volume with a link 
on the dual lattice, making world lines that define a conserved current density Jm- In 
Fig. |l](b) we show the three- volumes as squares since the time dimension is suppressed, and 
through the center of each square is the dual link associated with Jm- In order to isolate the 
solenoidal monopole currents, we construct the operator for the line integral of Jm around 
a dual plaquette, V x J M , from the four three- volumes (squares) shown in Fig. 0(b). Notice 
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from this construction that E and V x Jm take values at the same location within the unit 
cell of the lattice, indicated by the bold face line in Fig. [l|(b). 

IV. ANALYSIS 

In our previous work [0], we found that the confined phase of 4D U(l) lattice gauge 
theory exhibited a response to a static qq pair that could be described in terms of a dual 
version of the London theory ||. In this model, one combines the dual version of Ampere's 
law 

- cV x E = J M . (4) 

with the dual London equation governing the generation of the persistent currents in the 
monopole condensate 

VxJ M = ^E, (5) 

where A = (mc 2 /n s e 2 ) 1 / 2 is the London penetration depth, with e, m and n s the charge, mass 
and number density of the monopole condensate. The dual London theory also predicts, 
based on the single-valuedness of the condensate order parameter, the fluxoid quantization 
relation 

J E-dS-— j> J M - d£ = n$ e , (6) 

where n is an integer and $ e = e/Vhc is the quantum of electric flux. 

While this model explained our U(l) results, our results for SU(2) do not obey Eqn. 
Figure || shows the result of 480 measurements at /3=2.4 of E z and (V x Jm)z measured 
midway between the qq pair. Since -(V x Jm)z has a positive value at one point off axis, 
there is clearly no linear combination of the off axis data for E z and (V x Jm)z that will 
satisfy Eqn. (||). 

Therefore, we have adopted a dual form of Ginzburg-Landau (G-L) theory 0, which 
generalizes the London theory to allow the magnitude of the condensate density to vary in 
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space. As before, the superconducting order parameter is a complex function i[)(x), where 
|^(af)| 2 is the condensate density. We define ip(x) = y/n^f(x) exp(ia(x)), where n s is the 
London (bulk) condensate density, and / and a are real functions describing the spatial 
variation of the condensate; the London model presumes / = 1. The characteristic scale 
over which the condensate density varies is £, the G-L coherence length. 

The dual London Eqn. (y) is replaced with a dual version of the G-L equations || 

E = V x A E 

Ju = ^P{A E - |va) (7) 

and the condensate density obeys 

= -e 2 V 2 / + e 2 (Va - fA E ) 2 f -f + f, (8) 

while the fluxoid quantization relation of Eqn. (||) is generalized to 

J E-dS-^j j^J M ■ d£ = n$ e (9) 

The dual Abrikosov vortex for a flux tube along the z axis is a solution to Eqns. ([7p 
and (||) where the GL order parameter in polar coordinates varies as r/> = f(r) e ld with f(r) 
given approximately |Q by 

f(r) =tanh(0.9r/0- (10) 

To fit our data, we chose an analytic function for the azimuthal component of the 
monopole current of the form 



7 f \ ai 

J MKT) = — 
r 



;i + a 2 r)e" a2r - (1 + a 3 r)e- a3r (11) 



with a 1; a 2 , and a 3 as fitting parameters. The form of Eqn. ([□]) was simply chosen to 
ensure that Jm{t) vanished linearly at r = and faster than 1/r as r — > oo. The curl of the 
monopole current is then found from (V x Jm{t))z — {^-/i r )d(rJM)/dr and the parameters 
where fixed by fitting the data for the curl of the monopole current. The electric field data 
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was then fit using Eqns. ([U]) and (|TT|) in Eqn. (|7|). We excluded the point at r = in the fit 
to E z (r) since we are using a continuum approximation for the electric field to fit lattice data 
for the electric flux through a plaquette, and at the origin our continuum approximation to 
E z (r) diverges while the electric flux though the plaquette is still finite. 

We evaluate the fluxoid $ e by two distinct methods. The first identifies <3> e with $ tot , the 
net flux through the entire lattice. The second method comes from evaluating the fluxoid 
relation on the central plaquette (r = 0). If we interpret the data point E z (0) as actually 
representing the flux through a circle of unit area (radius 1/a/7t), Eqn. (??) gives 

A 2 2irrJ M (r 



E z (0) + 



(12) 



c f 2 (r) 

The curves in Figs. 0(a) and |2|(b) represent our fit to the data for (3=2.4. The resulting 
current distribution as a function of the distance from the qq axis is shown in Fig. |3|. 
The fitted values we obtained are A/a = 1.05 ± 0.12 and £/a = 1.35 ± 0.11. This latter 
result agrees with the naive estimate that £ is the distance where V x Jm changes sign, 
i.e., ~ 1.5a. The fluxoid found from Eqn. (|12|) was 0.16 ± 0.04, while that found from the 
total flux was §tot = 0.176 ± 0.003. The errors quoted are the statistical variance in the 
fitting parameters when the 480 measurements were divided into 4 sets of 120 points and 
fitting each set independently. For (3 = 2.5 we accumulated 208 measurements, which were 
grouped into 4 independent sets of 52 measurements for analysis. We found a penetration 
depth A/a = 1.59 ± 0.31, a coherence length £/a = 1.16 ± 0.08, and two flux quantum 
estimates as $ e = 0.18 ± 0.03 and $ to4 = 0.253 ± 0.005. 



V. CONCLUSIONS 

From this study we see that the dual superconductivity of the vacuum in SU(2) has a 
richer spatial structure than that seen in U(l), with a clearly apparent normal core region in 
the flux tube. The value of £ measures in effect the thickness of the boundary between the 
external, confining vacuum and the deconfined interior of a hadron. It will be interesting to 
investigate the structure of the flux tube in SU(3) lattice gauge theory. 
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The values we find for £ and A explain results obtained by Ivanenko et al. Jl3| on the 
existence of monopoles in SU(2) lattice gauge theory. They found the string tension between 
a qq pair correlated not with monopoles found by the deGrand-Toussaint construction ap- 
plied to elementary three-volumes, but rather to "extended" monopoles defined over several 
three- volumes. Since the correlation length £ measures the minimum length scale over which 
a well-defined monopole condensate exists, £ m a means the global monopole condensate is 
best measured over a larger region than an elementary three-volume. 

The gross spatial structure of the flux tube is fixed by the dimensionless Ginzburg-Landau 
parameter k = A/£. For a type II superconductor (k > l/v2) the flux tube is compact, 
while a type I superconductor (k < l/\/2) has a highly ramified flux structure [17]] in the 
form of corrugated cylinders or highly reticulated walls. Our result (re = 1.4 ± 0.2) is close 



to the borderline. We should point out that an analysis by Maedan et al. [IS] of an effective 
Lagrangian for confinement yields a value of k ^> 1 for SU(2). However, since our value for 
k depends on the size of the Wilson loop and on the scaling behavior of k(/3), we cannot yet 
reliably compare our results to a continuum model. Further analysis using larger Wilson 
loops will be needed to definitively establish the magnitudes of A and £. Rather, the goal 
of this work was to find an operator (V x Jm) that gave a clear signal of the nature of the 
confinement and showed the structure of the flux tube; an analysis of the electric flux alone 



does not elucidate this structure [19 
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FIGURES 

FIG. 1. Operators for (a) electric field and (b) V x Jm on a fixed time slice. 

FIG. 2. Profile of (a) E z and (b) V x Jm as a function of distance from the qq axis for f3 = 2.4, 
together with the best fit. 

FIG. 3. Profile of the current distribution Jm derived from the fit to the data of Fig. |2[ 
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